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Abstract 

The relations between smooth and peaked soliton solutions are reviewed for the Camassa- 

Holm (CH) shallow water wave equation in one spatial dimension. The canonical Hamiltonian 
formulation of the CH equation in action-angle variables is expressed for solitons by using 
the scattering data for its associated isospectral eigenvalue problem, rephrased as a Riemann- 
Hilbert problem. The momentum map from the action-angle scattering variables T*{T^) to 
the flow momentum (X*) provides the Eulerian representation of the A^-soliton solution of 
CH in terms of the scattering data and squared cigcnfunctions of its isospectral eigenvalue 
problem. The dispersionless limit of the CH equation and its resulting peakon solutions 
are examined by using an asymptotic expansion in the dispersion parameter. The peakon 
solutions of the dispersionless CH equation in one dimension are shown to generalize in 
higher dimensions to peakon wave-front solutions of the EPDiff equation whose associated 
momentum is supported on smoothly embedded subspaces. The Eulerian representations of 
the singular solutions of both CH and EPDiff are given by the (cotangent-lift) momentum 
maps arising from the left action of the diffeomorphisms on smoothly embedded subspaces. 
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1 Shallow water background for the CH equation 

Euler's equations for irrotational incompressible ideal fluid motion under gravity with a free surface 
have an asymptotic expansion for shallow water waves that contains two small parameters, e and 
S^, with ordering e > S^. These small parameters are e = a /ho (the ratio of wave amplitude 
to mean depth) and 5^ = {ho/lxY (the squared ratio of mean depth to horizontal length, or 
wavelength). Euler's equations are made non-dimensional by introducing x = l^x' for horizontal 
position, z = h^z' for vertical position, t = {lx/co)t' for time, rj = at]' for surface elevation and 
If = {glxa/ Cojif' for velocity potential, where cq = V 9^0 is the mean wave speed and g is the 
constant gravity. The quantity a = a' / {hopcD is the dimensionless Bond number, in which p is 
the mass density of the fluid and a' is its surface tension, both of which are taken to be constants. 
After dropping primes, this asymptotic expansion yields the nondimensional Korteweg-de Vries 
(KdV) equation for the horizontal velocity variable u — (px{x,t) at linear order in the small 
dimensionless ratios e and 5^, as the left hand side of 

3e 52 

Ut + Ux + —uux + —(1 - ^a)uxxx = 0(e5^) . (1) 

Here, partial derivatives are denoted using subscripts, and boundary conditions are u — Q and 
Ux = ^ bX spatial infinity on the real line. The famous sech^{x — ct) traveling wave solutions 
(the solitons) for KdV (1) arise in a balance between its (weakly) nonlinear steepening and its 
third-order linear dispersion, when the quadratic terms in e and 5^ on its right hand side are 
neglected. 
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On the right hand side of equation (1), a normal form transformation due to Kodama [40] has 
been used to remove the other possible quadratic terms of order O(e^) and 0{S^). The remaining 
quadratic correction terms in the KdV equation (1) may be collected at order 0{e5'^). These terms 
may be expressed, after introducing a "momentum variable," 

m — u — uS^Uxx , (2) 

and neglecting terms of cubic order in e and 6^, as 

e 5^ 
mt + mx + -{unij: + bmux) + —(1 - Sa)uxxx = . (3) 

In the momentum variable m — u — v^'^Uxxi the parameter v is given by [21] 

19- 30(7 -45(72 
60(1-3(7) ■ 

Thus, the effects of 5^- dispersion also enter the nonlinear terms. After restoring dimensions in 
equation (3) and rescaling velocity m by (6+1), the following "6-equation" emerges, 

mt + c^vfix + umx + h mu^ + Vu^xx = , (5) 

where m — u — o?Uxx is the dimensional momentum variable, and the constants o? and F/cq are 
squares of length scales. When o? — > 0, one recovers KdV from the 6-equation (5), up to a rescahng 
of velocity. Any value of the parameter h ^ —\ may be achieved in equation (5) by an appropriate 
Kodama transformation [21]. 

As we have emphasized, the values of the coefficients in the asymptotic analysis of shallow water 
waves at quadratic order in their two small parameters only hold, modulo the Kodama normal- 
form transformations. Hence, these transformations may be used to advance the analysis and 
thereby gain insight, by optimizing the choices of these coefficients. The freedom introduced by the 
Kodama transformations among asymptotically equivalent equations at quadratic order in e and 
also helps to answer the perennial question, "Why are integrable equations so ubiquitous when 
one uses asymptotics in modeling?" Namely, there may be special values of the free parameters in 
the normal-form transformations of the asymptotics that allow integrability. 



Integrable cases of the 6-equation (5). The cases 6 = 2 and 6 = 3 are special values. For 
these values, the 6-cquation becomes completely integrable as a Hamiltonian system. For 6 = 2, 
equation (5) specializes to the integrable CH equation of Camassa and Holm [9]. The case 6 = 3 in 
(5) recovers the integrable DP equation of Degasperis and Procesi [19]. These two cases exhaust 
the integrable candidates for (5), as was shown using Painleve analysis. The 6— family of equations 
(5) was also shown in [46] to admit the symmetry conditions necessary for integrability, only in 
the cases 6 = 2 for CH and 6 = 3 for DP. 

The 6-equation (5) with 6 = 2 was first derived in Camassa and Holm [9] by using asymptotic 
expansions directly in the Hamiltonian for Euler's equations governing inviscid incompressible flow 
in the shallow water regime. In this analysis, the CH equation was shown to be bi-Hamiltonian 
and thereby was found to be completely integrable by the inverse scattering transform (1ST) on 
the real line. This development of 1ST for CH equation (5) with 6 = 2 is discussed in Section 2. 

Camassa and Holm [9] also discovered the remarkable peaked soliton (peakon) solutions of 
(58,59) for the CH equation on the real line, given by (5) in the case 6 = 2. The peakons arise as 
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solutions of (5), when cq — and F = in the absence of hnear dispersion. Peakons move at a 
speed equal to their maximum height, at which they have a sharp peak (jump in derivative). The 
single peakon solution is 



Unlike the KdV soliton, the peakon speed is independent of its width (a). Periodic peakon solutions 
of CH were treated in Alber et al. [3]. There, the sharp peaks of periodic peakons were associated 
with billiards reflecting at the boundary of an elliptical domain. These billiard solutions for the 
periodic peakons arise from geodesic motion on a tri-axial ellipsoid, in the limit that one of its 
axes shrinks to zero length. 

The CH equation was found after its derivation as a shallow water equation in [9] to fit into 
a class of integrable equations derived previously by using hereditary symmetries in Fokas and 
Fuchssteiner [23] . See Fuchssteiner [25] for an insightful history of how the shallow water equation 
(5) in the integrable case with h — 2 relates to the mathematical theory of hereditary symmetries. 

Equation (5) with 6 = 2 was recently re-derived as a shallow water equation by using asymptotic 
methods in three different approaches in Dullin et al. [21], in Fokas and Liu [24] and also in Johnson 
[37]. These three derivations all used different variants of the method of asymptotic expansions 
for shallow water waves in the absence of surface tension. Only the derivation in Dullin et al. 
[21] took advantage of the parametric freedom in the Kodama normal-form transformations of the 
asymptotic expansion results at quadratic order. 

The effects of the parameter h on the solutions of equation (5) were investigated in Holm and 
Staley [35], where h was treated as a bifurcation parameter, in the limiting case when the linear 
dispersion coefficients are set to cq = and F = 0. This limiting case allows a variety of special 
solutions for different ranges of the values of 6, in which the two nonlinear terms in equation (5) 
balance each other in the absence of linear dispersion. 

Since its first appearance in [9] , the CH equation has been the centre of a confluence of scientiflc 
endeavors that includes water waves, integrable systems, PDE analysis, asymptotics, geometry and 
Lie groups. This confluence has led to a continued interest and opportunities for contributions 
from many different fields in mathematics. The interest in CH solutions may be measured by 
noticing that the entry 'peakon' has acquired many thousands of Google hits. 

Plan of the paper 

This paper aims to review some of the geometric highlights of recent work on the CH equation. 
It is certainly not exhaustive. It mainly focuses on comparing the soliton theory for smooth CH 
solutions with the peakon theory for its singular solutions in the dispersionless limit. Section 2 
briefiy explains the application of the method of the inverse scattering transform (1ST) for obtaining 
the soliton solutions of CH equation. The set of scattering data is introduced and the formulation 
of the inverse scattering as a Riemann-Hilbert problem is outlined. The solution is expressed via 
the scattering data in a form that admits the peakon limit in the sections that follow. In Section 3 
the map between the action-angle variables (expressed via the scattering data) and the momentum 
of the CH solution is formulated as a momentum map from the symplcctic action-angle variables 
to the dual of the Lie algebra of smooth vector fields on the real line. This is a Poisson map, but 
the noncanonical bi-Hamiltonian structure of the CH equation which led to the discovery of its 
isospectral problem in [9] and its geometrical significance as geodesic flow on the diffeomorphism 
group are not discussed here. In Section 4 we introduce the peakons as singular solutions that 
appear in dispersionless limit. The A^-peakon solution is governed by a flnite-dimensional integrable 




(6) 
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dynamical system. Section 5 presents the multi-peakon solution as a limiting case of the CH 
multi-soliton solution and points out the similarity between the dynamics of the peakon system 
and the well known Toda lattice. In Section 6 we comment briefly on the existence of additional 
integrals of motion of the peakon system, a property known as superintegrability. The two-peakon 
system is analysed explicitly. Section 7 presents the compacton and pulson solutions, which are 
singular solutions, similar to peakons that can also be represented by the singular momentum 
map for peakons discussed in [33]. In Section 8 we deal with generalizations of the CH equation 
in higher dimensions. These also allow for singular solutions. Although such multidimensional 
generalizations are unlikely to be integrable, numerical studies show that their solutions are stable 
and interact elastically in collisions. In Section 9 we comment briefly on some limitations of our 
present discussion and point out three open problems for further research. 

2 Soliton Solutions of CH equation from Inverse Scatter- 
ing 

2.1 Inverse scattering for the KdV equation 

One of the most significant results in the theory of nonlinear partial differential equations was the 
development of Gardner, Greene, Kruskal and Miura [28, 29] of a method for the exact solution 
of the initial-value problem for the KdV equation. Prior to their work the only known exact 
solutions of KdV were the travelling wave solutions. The idea is based on a representation of the 
equation {ut + Guu^ + u^xx = 0) as a compatibility condition of two eigenvalue problems with a 
time-independent spectral parameter A: 



where 7 is an arbitrary constant. The method is conceptually analogous in many ways to the 
Fourier transform method for solving linear equations. This method is now known as Inverse 
Scattering Transform (1ST). It was recast by Lax [41] in a general framework that allows the 
1ST to be used for solving other nonlinear PDEs. An important consequence is the relationship 
between the discrete eigenvalues and the characteristics of the solitons that emerge asymptotically. 
Another feature of soliton collisions is the preservation of soliton identities after the interaction 
(asymptotically in time). For example, collisions of KdV solitons only result in a phase shift from 
the positions where they would have been without the interaction. 

The KdV was formulated as a completely integrable Hamiltonian system in a work by Faddeev 
and Zakharov [55], the Hamiltonian form was also noted by Gardner [27]. 

Reviews of 1ST may be found, for example, in Ablowitz et al. [1], Dubrovin et al. [20], Novikov 
et al. [47]. For discussions of other related bi-Hamiltonian equations, see [19]. 

2.2 Inverse scattering for CH solitons with dispersion 

In this section we outline the application of the 1ST for the CH equation. We use the form with a 
linear dispersion term 



*t = (ii^ + 7)*- (4A + 2w)^'^ 



(7) 
(8) 



'^t ~ Uxxt + 'iuJUx + 3uUx — 




(9) 
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where a; is a real constant. The equation in this form also appears as a model of unidirectional 
propagation of shallow water waves over a flat bottom [9, 37, 38] as well as that of axially symmetric 
waves in a hyperelastic rod [18]. It can be obtained from the 6-equation (5) with b — 2 via Galilean 
transformation that removes the u^^;,. term. 

Equation (9) admits a Lax pair formulation [9] 

= J(l-(m + a;))^' (10) 

% = -0 + «^*. + y^ + 7^ (11) 

where 7 is an arbitrary constant. We will use this freedom for a proper normalization of the 
eigenf unctions. 




Figure 1: A Gaussian initial condition for the CH equation breaks up into an ordered train of 
solitons as time evolves (the time direction being vertical). The soliton train eventually wraps 
around the periodic domain, thereby allowing the leading solitons to overtake the slower emergent 
solitons from behind in collisions that cause phase shifts, as discussed in [9]. 

In our further considerations m will be a Schwartz class function, u > and m{x, 0) + u > 0. 
Then m{x,t) + u > for all t [11]. Let us introduce a new spectral parameter k, such that 

x{k)^^(i+ik^y (12) 

The spectrum of the problem (10) under these conditions is described in [10]. The continuous 
spectrum in terms of k corresponds to A; - real. The discrete spectrum (in the upper half plane) 
consists of finitely many points kn — iHn, n = 1, . . . , N where Kn is real and < k„ < 1/2. 

For all real A; 7^ a basis in the space of solutions of (10) can be introduced, ip{x, k) and -0(0;, k), 
fixed by its asymptotic when x ^ 00 [10]: 

'^{x.k) = e-*^^ + o(l), x^oo (13) 

Another basis can be introduced, ip{x, k) and (p{x, k) fixed by its asymptotic when x — >■ —00: 

^{x,k) = 6-^^=^ + 0(1), x^-00 (14) 
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and the relation between the two bases is [13] 

ip{x, k) = a{k)ip{x, k) + b{k)ip{x, k). (15) 

where 

|a(A;)|'- |6(A;)|' = 1. (16) 

The quantity lZ{k) = h{k)/a{k) is called the reflection coefficient, without any danger of con- 
fusing h{k) with the bifurcation parameter h introduced in (5). 

All of the required information about the scattering, i.e. a{k) and 6(/c), is provided by IZik) for 
A; > only [15]. It is sufficient to know TZ{k) only on the half hue A; > 0, since a{k) = a{—k), 
l{k) = h{-k) and thus n{-k) = n{k). 

The constant 7 in (11) can be chosen for each eigenfunction in such a way that a{k) does not 
depend on t and is a generating function of the integrals of motion [15]. At the points of the 
discrete spectrum, a{k) has simple zeroes [10], therefore tp and are linearly dependent (15): 

</7(x, IKn) = hnil{x, -iK„). (17) 

In other words, the discrete spectrum is simple, there is only one (real) eigenfunction 
corresponding to each eigenvalue iKm and we can take this eigenfunction to be 

^^''\x) = ^{x,iKn) (18) 

The asymptotic behaviour of (/?^"\ according to (14), (17) is 

^i^){x) = e""^ + o(e"""), x^-oo; (19) 
<^(")(x) = 6„e-''"^ + o(e-''"^), x^oo. (20) 

The sign of 6„ obviously depends on the number of the zeroes of (p^'^\ Suppose that 

< Ki < K2 < . . . < kn < 1/2. (21) 

Then from the oscillation theorem for the Sturm-Liouville problem [7], (p^^^ has exactly n — 1 
zeroes. Therefore 

bn^{-ir-'\bn\. (22) 

The set 

S = {n{k) {k>0), Kn, Rn = bJia'{tKn), u = 1, . . . N} (23) 

is called the scattering data. The Hamiltonians for the CH equation in terms of the scattering data 
are presented in [15]. The time evolution of the scattering data is [13]: 

TZ{k,t) = 7^(A;,0)e-'?^ (24) 
Kit) = i?„(0)exp(^^t^, (25) 

where Rn{t) is always a positive quantity [13] and A„ = A(fi;„). 



D.D. Holm and R.I. Ivanov 



Smooth and Peaked CH Solitons 



8 



The scattering coefficient a{k) is analytic for Im A; > with asymptotic [15] 

e'f^''a{k) 1, \k\ oo. (26) 
where ^ is a constant (integral of motion): 



The asymptotics of the Jost solutions (or, rather the following quantities, depending on the 
Jost solutions) for \k\ oo have the form [13] 

l{x,k) = V'(x,A;)e^=Xo(x) + ^ + ^ + ..., (28) 

^{x,k) ^ ^{x,k)e''^^^^^^Xo{x) + ^ + ^ + ... (29) 
— k k'^ 



where Xq{x) = [uj/{m{x) + u})]^^"^ > 0, 



y{x) = K 



x+ {Xo^{x) - l)dx 

J oo 



(30) 



Moreover, (28) is analytic for Im A; < 0, (29) is analytic for Im A; > [13]. 
Prom (15) with (28) and (29) we obtain 



e^^!^a{k) 



ilj{x, k) + 7^(A;)^(a;, A;)e'*'=^w/V'^. (31) 



The function v?(x, k) / [e^^^ a{k)) is analytic for Im A; > 0, ^(x, k) is analytic for Im A; < 0. Thus, 
(31) represents an additive Riemann-Hilbert Problem (RHP) with a jump on the real fine, 

given by n{k)ip{x,k)e^'^y^''^/^. 

The solution of this RHP follows a standard pattern, see [13] for more details. The solution 
simplifies considerably in the case TZ{k) = 0. This is the iV-soliton solution: 

RJt)e-^^^y/^'i/j(x, -iKn) 
ij(x,k) = Xo(x)-tY^^ : =^ -, lmk<0, 32 

OfVfl IS, 

n=l 

= X,{x)-iY,— . Imfc>0. (33) 

From (32) one has a linear system for the quantities t/'(a;, —inmt) with a solution 

N 

±{x,-iKn,t)^XQ{x)Y,Kl[yA, n^l,...,N. (34) 

p=i 

where 
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n,p=l '^^ 2 



The substitution «; = i/2 in (33) with a (i/2) = 1 gives 



From (36) and (37) there follows a parametric representation 



X = X{y,t) = ^+\n^, 

J- 



h - 1-E -^--'b.'i- 

2.3 Parametric form of the dispersive CH soliton solution 

From (30) and (38) one can compute the solution in parametric form 

u{X{y,t),t)^Xt{y,t), x^X{y,t), 

where X{y,t) is given in terms of the scattering data in (38), (39). 
Upon introducing new notations 

^ f y 2t 
= 2.,^-— + -+.,0 

= ^ln5M 

^, = In 1-2"'- 



1 + 2k 

la = In 



■J 

2 



ACj ~l~ ACj 

one can rewrite the expression for j± (39) in the form [44, 48, 49, 50] 
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N 

1=1 l<i<j<j 



f± = 

A»=0,1 

The solution for m can be obtained from (40). First we notice that 

du{X{y,t),t) 



(45) 



dy 

and also 



u^{X{y,t),t)Xy. 



du{X{y,t),t) ^ dX, 
dy dy 

Thus Ua:{X{y,t),t) = Xty/Xy. Similarly, 

«..(X(|/,t),t) = i-(^^) 



and 



1 /X, 



m{X{t,y),t) = u{X{t,y),t)-u,,{X{t,y),t)^Xt-^ ( 

/oo 
P{y,t)5{x-X{y,t))dy 
-oo 



(46) 



with p(^,t)=x,X,-(^] , (47) 

y / y 



1 /"^ 

u{x,t) = -y P(|/,t)exp(-|a;-X(|/,t)|)dy. (48) 

As we shall see in Section 8, this representation is also useful in the study of multidimensional 
solutions. 

2.4 Relation to KdV hierarchy 

The spectral problem (10) is gauge equivalent to a standard Sturm-Liouville problem, well known 
from the KdV hierarchy, cf. (7), with short notation q = m + u: 

A 1 

-%y + U{y)^ = A^=4-^' 

m = q'^'^, %-V~q. (49) 

45(y)^4?(y) IQq^y) ^ ^ 

Note that (49) leads to two possible expressions for the change of the variables in the Liouville 
transformation: 

/X 
{■\/q{x') — ^/u)dix' + const, (51) 
-oo 

l-X 

y — ^/ux-\- I (-\/g(x^ — -\/^)dx' + const. (52) 

J oo 
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These two possibihties, (51), (52) are only consistent iff 

/oo 
iV q{x) ~ Vuj)dx — constant, 
-oo 

which is always the case, since the integral under question is (up to a multiplier) the Casimir 
function P (27). 

The matching of the CH hierarchy to KdV hierarchy requires solving the Ermakov-Pinney 
equation (50) [10] , which is not straightforward and leads to the same solution (40) in parametric 
form. 



3 Momentum map formulation with action-angle variables 

The canonical Poisson brackets for the scattering data of the CH equation are computed in [15] 
where also the action-angle variables are expressed in terms of the scattering data. Let us consider 
the action variable for the N-soliton solution. (These considerations can be extended easily to the 
variables of the continuous spectrum.) The angle variable is = lni?„(t); it is linear in time t 
and $„ = 4:K.n/Xn- Let us introduce A„ := Ann/Xn into Hamilton's principle 5S = 0, with 

%,$„,n,]= / U[u] + Y.^n[K- K[u]))dt. 

^ n=l ^ 

Here the Lagrange multipher H„ enforces the action-angle relation for the CH scattering data as a 
shift of the angles at constant angular frequencies A„, with n — 1, . . . ,N . The stationary 
variation is 

N N \ 

dt. 

n=l n=l 



o=ss= I u^-j2^j-^yu+j2{^n-AMyiin-j2^j^- 

\ ^ r> — 1 ^ n — 1 n — 1 



Since by definition m — 5i/5u is the momentum, 5S — implies the Eulerian representation 

n=l 

^n^K[u\ and rin = 0. (54) 
Relation (53) is the momentum map 

($„, n„) G T*T^ -^meX* 

for the toral action (54) on the angles at constant angular frequencies A„. This momentum 
map from the action-angle scattering variables T*(T^) to the flow momentum X*(M) (dual to 
the smooth vector fields X(]R) on the real line) provides the Eulerian representation of the A^- 
soliton solution of CH in terms of the scattering data and squared eigenfunctions of its isospectral 
eigenvalue problem. Momentum maps for Hamiltonian dynamics are reviewed in [43], for example. 

By using the spectral quantities of the A^-soliton solution (recall: A„ = (1 — 4K,f^)/u) one may 
express the variation of the spectrum with respect to the CH solution in terms of the squared- 
eigenfunctions of the isospectral problem as [15] 

^"^^ - ^^ + ^'''^-i?n(^)[^(x,-^«„,^)]^ 



6m{x,t) 20JKn\n 
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in which ip{x, —iKn,t) is the eigenfunction that belongs to eigenvalue A^, see (17). 

On the other hand, the expansion of u{x, t) over squares of eigenfunctions is given by [14] 



Consequently, 



or 



n=l 



^ 4k. — 

m{x,t) = J2 u^^n{t){i - d^)[i^{x, -m„,^)]^ 

n=l 



N 

m{x,t) = ^n„J„(x,i) , 



n=l 



where n„ and Jn{x,t) denote explicitly 



= ,''' Rnmi-d'Mx,-ZKn,t)]\ (57) 

Thus, the momentum map (53) from the action-angle variables under going dynamics (54) to 
the Eulerian representation of the momentum for the CH solution is expressed in terms of the 
scattering data and squared eigenfunctions of its A'"-soliton isospectral eigenvalue problem. Perhaps 
not unexpectedly, this momentum map may be applied to the action-angle representation of the 
solution of any integrable Hamiltonian PDE. 



4 Peakons 

4.1 Peakons: the singular solution ansatz 

Camassa and Holm [9] discovered the "peakon" solitary travehng wave solution for a shallow water 
wave, 

M(x,t) = ce~l"-'^*l/", (58) 

whose fluid velocity u is a function of position x on the real line and time t. The peakon traveling 
wave moves at a speed equal to its maximum height, at which it has a sharp peak (jump in 
derivative). Peakons are an emergent phenomenon, solving the initial value problem for a partial 
differential equation derived by an asymptotic expansion of Euler's equations using the small 
parameters of shallow water dynamics. Peakons are nonanalytic solitons, which superpose as 



^(^' t) = I \ EPa(i)^(x - qa{t))/a) , 

a=l a=l 



(59) 
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for sets {p} and {q} satisfying canonical Hamiltonian dynamics. Peakons arise for shallow water 
waves in the Umit of zero linear dispersion in one dimension. Peakons satisfy a partial differential 
equation (PDE) arising from Hamilton's principle for geodesic motion on the smooth invertible 
maps (diffcomorphisms) with respect to the Sobolev norm of the fluid velocity. Peakons 
generalize to higher dimensions, as well. We explain how peakons were derived in the context of 
shallow water asymptotics and describe some of their remarkable mathematical properties. 

Peakons were first found as singular soliton solutions of the completely integrable CH equation. 
This is equation (5) with b — 2, now rewritten in terms of the velocity, as 

Ut + CqUx + Suux + ruxxx 

= a^{u^xt + '^u^u^^ + uu^xx) ■ (60) 

Peakons were found in [9] to arise in the absence of linear dispersion. That is, they arise when cq = 
and r = in CH (60). Specifically, peakons are the individual terms in the peaked A^— soliton 
solution of CH (60) for its velocity, in the absence of hnear dispersion. Each term in the sum (59) 
is a solition with a sharp peak at its maximum. Hence, the name "peakon." Expressed using its 
momentum, m = (1 — a'^d'^)u, the peakon velocity solution (59) of dispersionless CH becomes a 
sum over a delta functions, supported on a set of points moving on the real line. Namely, the 
peakon velocity solution (59) implies 

N 

m{x, t) = Pait)5{x - g„(t)) , (61) 

o=l 

because of the relation (1 — Q;^9^)e~l^l/" = 2aS{x). These solutions satisfy the 6-equation (5) for 
any value of b, provided Cq — and F = 0. As we shall discuss later, the peakon momentum 
relation (61) is again a momentum map. 

Thus, peakons are singular momentum solutions of the dispersionless 6-equation, although they 
are not stable for every value of b. From numerical simulations [35], peakons are conjectured to be 
stable for 6 > 1. In the integrable cases 6 = 2 for CH and 6 = 3 for DP, peakons are stable singular 
soliton solutions. The spatial velocity profile e~l^l/"/(2Q;) of each separate peakon in (59) is the 
Green's function for the Helmholtz operator on the real line, with vanishing boundary conditions 
at spatial infinity. Unlike the KdV soliton, whose speed and width are related, the width of the 
peakon profile is set by its Green's function, independently of its speed. 

4.2 Integrable peakon dynamics of CH 

Substituting the peakon solution ansatz (59) and (61) into the dispersionless CH equation, 

rrit + wrix + 2mUx — , with m — u — Q?Uxx , (62) 

yields Hamilton's canonical equations for the dynamics of the discrete set of peakon parameters 
Pa{t) and qa{t), 

for a = 1, 2, . . . , A?", with Hamiltonian given by [9], 

1 ^ 

a,b=l 
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Or explicitly, 

qa = -5]p,e-l^»-^^l/" (65) 

6=1 

N 

= f^5^P6e-l''"-^''l/"sgn(g„-gfe) (66) 

^ 6=1 

Thus, one finds that the points x = qa{t) in the peakon solution (59) move with the fiow of 
the fluid velocity u at those points, since u{qa{t),t) = qa{t). This means the qa{t) are Lagrangian 
coordinates. Moreover, the singular momentum solution (61) is the Lagrange-to-Euler map for 
an invariant manifold of the dispersionless CH equation (62). On this finite-dimensional invariant 
manifold for the partial differential equation (62), the dynamics is canonically Hamiltonian. 

With Hamiltonian (64), the canonical equations (63) for the 2N canonically conjugate 
peakon parameters Pa{t) and qa{t) were interpreted in [9] as describing geodesic motion on the 
A?"— dimensional Riemannian manifold whose co-metric is g^^{{q}) — e"'^""^''!/". Moreover, the 
canonical geodesic equations arising from Hamiltonian (64) comprise an integrable system for any 
number of peakons N. This integrable system was studied in [9] for solutions on the real line, and 
in [3, 16] and references therein, for spatially periodic solutions. 

The integrals generated by the action variables in terms of the coordinates can be recovered as 
tr(L'^), where L is the Lax operator for the peakon system [9]. 

Being a completely integrable Hamiltonian soliton equation, the continuum CH equation (60) 
has an associated isospectral eigenvalue problem, discovered in [9] for any values of its dispersion 
parameters Cq and F. Remarkably, when Cq = and F = 0, this isospectral eigenvalue problem has 
a purely discrete spectrum. Moreover, in this case, each discrete eigenvalue corresponds precisely to 
the time-asymptotic velocity of a peakon. This discreteness of the CH isospectrum in the absence 
of linear dispersion implies that only the singular peakon solutions (61) emerge asymptotically in 
time, in the solution of the initial value problem for the dispersionless CH equation (62). This is 
borne out in numerical simulations of the dispersionless CH equation (62), starting from a smooth 
initial distribution of velocity [26, 35]. 

Figure 2 shows the emergence of peakons from an initially Gaussian velocity distribution and 
their subsequent elastic collisions in a periodic one-dimensional domain. This figure demonstrates 
that singular solutions dominate the initial value problem and, thus, that it is imperative to go 
beyond smooth solutions for the CH equation; the situation is similar for the EPDiff equation. 

Peeikons as mechanical systems. Governed by canonical Hamiltonian equations, each 
A?"— peakon solution can be associated with a mechanical system of moving particles. Calogero 
et al. [8] further extended the class of mechanical systems of this type. The r-matrix approach 
was applied to the Lax pair formulation of the A^— peakon system for CH by Ragnisco and Br- 
uschi [52], who also pointed out the connection of this system with the classical Toda lattice. A 
discrete version of the Adler-Kostant-Symes factorization method was used by Suris [53] to study 
a discretization of the peakon lattice, realized as a discrete integrable system on a certain Poisson 
submanifold of gl{N) equipped with an r-matrix Poisson bracket. Beals et al. [6] used the Stieltjes 
theorem on continued fractions and the classical moment problem for studying multi-peakon so- 
lutions of the CH equation. Generalized peakon systems are described for any simple Lie algebra 
by Alber et al. [3]. 
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Figure 2: A Gaussian initial condition for the CH equation breaks up into an ordered train of 
peakons as time evolves (the time direction being vertical). The peakon train eventually wraps 
around the periodic domain, thereby allowing the leading peakons to overtake the slower emergent 
peakons from behind in collisions that cause phase shifts, as discussed in [9]. 



5 Peakon limit of the CH soliton solutions 



The hmit a; — > in the A^-soliton solution u{x, t) produces the A^^-peakon solution (59). The limiting 
procedure is described in detail in [45]. Due to (12) one can write for the discrete eigenvalues 



2Kj = = 1 



1 



-(jjXj + . . . . 



(67) 



The solution (38) depends explicitly on kj (67) and the limit can be computed with (67) by 
taking u ^ and keeping the eigenvalue Xj constant. The result is the expression (59) with 



Pi 



4^(0) ^(2) 
^(1) ^(1) ■ 



Qi — aln 



2D 



(0) 



D 



(2) 
N-i 



(i = l,2,...,Ar), 
(i = l,2,...,Ar), 



(68) 
(69) 



where 



l<ii<i2<-:<in<N 



n^l,2,...,N, 



A^{z,,Z2,...,Zn) = n (K-KJ', (n>2), (70) 



l<l<m<n 

2 . 



Ri(t) = Ri(0)e^i\ (xio = lni?,(0)). (71) 
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The quantities Dn"*^ are called Hankel determinants. By definition — 1. In general, 

Hankel determinant is a determinant of n x n matrix of the form = det(a^J^^) where Ojj are 
elements of a sequence, i.e. 

In this particular case 



N 



A, = ^Aii?,(t). (72) 



i=l 



Similcirity of peakon lattice and Toda lattice Hankel determinants appear in the solutions 
of other integrable systems, e.g. Toda lattice, e.g. see [31]. The Toda equation [54] 

dt ~ 

— = Pn, neZ 

is one of the most important integrable systems. We have a finite chain with N nodes under the 
fixed ends boundary conditions go = ~Qn+i — oo. The solution is in the form 

^(0) 

Quit) = gi(0) + ln-^, (73) 

qi{0) = --^lnDSJ^(O) is a constant. (74) 
The Hankel determinants are obtained by a similar sequence, 

N 



Ai = $^Aii?,(t), (75) 



1=1 



Ri{t) = i?i(0)e-^**, (76) 

where Aj are N different constants (eigenvalues of the Lax matrix) and the quantities i?i(0), 
i — 1, . . . , N represent another set N of constants. 

Due to its simple form, the A'"-peakon solution can be used as an approximation of the A^-soliton 
CH solution when the dispersion term is small (and the term 2ujUx can be neglected). Similarly, the 
Toda chain with complex dynamical variables (the so-called Complex Toda Chain - CTC) provides 
an approximation for the N-soliton solution of the Nonlinear Schroedinger Equation (NLS) 

1 , ,2 

lUt + -u^x +\u\ u = 0, 

see [32, 30, 31] for more details. Such an approximation is called adiabatic approximation and 
means that the A^-soliton solution consists of N well separated solitons 

2i/fce*P'"'=(^~^'=(*)+'''^(*)] 

u{x,t) ^ 2^ 



k-- 



^ cosh(2z/fc(a; - ik{t)) ' 
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i.e. the overlap between the sohtons is small. The variables qn{t) of the CTC are related to the 
NLS sohtons parameters by 

qk{t) = -2i'Qik{t) + i {2iiQikit) - Skit)) + const 

where ^fc, 5^, = |Cfc and Vk = (|^fc — nl)^^'^ characterize the center-of-mass position, the phase, 
velocity and amplitude respectively of the k-th soliton in the chain: z/q and are the average 
amplitude and velocity of the soliton train. and 6k can be obtained as the real and imaginary 
parts of qk{t). Such soliton trains and their asymptotic behavior appear to be important for the 
needs of soliton based fiber optics communications. 

6 Superintegrability of the peakon system 

Suppose we have an integrable system with 2N- dimensional phase space which in terms of the 
Action-Angle (canonical) variables can be represented as 

A„ = 0, 6„ = A„, n^l,2,...,N, (77) 

or, if there exists a bracket such that ^ 

AJ = ^6ni, $J = {A„ Ai} = 0, (78) 

the system is Hamiltonian with a Hamiltonian 

hN = Al + ...A%. (79) 

The integrals A„, n = 1,2, ... ,N are clearly in involution, which guarantees the integrability of 
the system. There is however another set of integrals: 

/. = ($ . - $,.+i)(Ai + . . . + Ajv) - (A,- - A,-+i)($i + . . . $iv) 

j = l,2,...,Ar-i. (80) 

If the 'action' variables A„ are all different, the set (80) is functionally independent from the set 
A„, n — 1,2, . . . , N. In addition, the integrals (80) form another set of A^-integrals in involution 
together with H. Due to the existence of two sets of functionally independent integrals in involution 
such systems are termed superintegrable. 

An example of such system is Toda lattice, see the discussion in [2]. The peakon system is also 
superintegrable. The canonical variables in terms of the scattering data for CH equation can be 
used in the peakon limit: A„ = 2/A„, = In Rn{t), the Hamiltonian (79) is also a peakon limit 
{u! — > 0) of the A^-soliton Hamiltonian [15] 

rr ( N 2V"/, 1-2k„ 4/t„(l + 4/t^) ^^ 

H.i.)^. X:(;in^-^ + ^^-^j. 

^The 1/2 coefficient in the definition of the bracket appears in order to match it to the Poisson bracket used in 
the CH Hamiltonian formulation. 
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The Poisson bracket for the CH peakon solution is 

/oo r /t A R 

— (md + dm)—dx. (81) 
dm dm 

and the scattering data satisfy (78) with respect to (81), see [15] for the details. 

More interesting are the integrals (80). From (69) and (68) one can recover these integrals in 
coordinate form. For example, when iV = 2 we have 



h = In ^J_ + P^-P- +^^Ml+^+i^Pl\ (82) 



\fj-pi+P2 P1+P2V Ot P2, 

where J = (pi - P2f + 4pip2e"l^^"^^l/°^. (83) 

Note that Ji depends on both combinations qi + q2 and qi — q2 as well as the momentum 
variables. The Hamiltonian /12 (which depends only on qi — q2 and the momentum variables) and 
Ii form a complete system of integrals in involution. The integration of the 2-peakon system with 
these integrals can be performed as follows. First, one can express qi and q2 in terms of Ji, /12 and 
the momentum variables: qi — qi(pi,P2, h, ho). Next, 

% . , % . /o.X 

(li = j^Pi + 75—^2- 84) 
opi dp2 

The substitution of qi, pi from (65) (66) to (84) produces an algebraic equation that gives, say p2 as 
a function of pi. Then (66) is an ODE for pi of the form pi = /(pi, /i, /?.2). Clearly, from practical 
viewpoint it is much more convenient to work with the other system of integrals in involution: /i2 
and the conserved momentum P — pi + P2- Note that J = 4/i2 — 3P^ is itself an integral. 



7 Other singular solutions: the dispersionless 6-equation 

Pulsons: Generalizing the peakon solutions of the dispersionless 6-equation for other 
Green's functions. The Hamiltonian in equation (64) depends on the Green's function for 
the relation between velocity u and momentum m. However, the singular momentum solution 
ansatz (61) is independent of this Green's function. Thus, as discovered in Pringer and Holm [26], 
The singular momentum solution ansatz (61) for the dispersionless equation, 

mt + umx + 2mUx = 0, with u = g *m, (85) 

provides an invariant manifold on which canonical Hamiltonian dynamics occurs, for any choice 
of the Green's function g relating velocity u and momentum m by the convolution u = g * m. 

The fluid velocity solutions corresponding to the singular momentum ansatz (61) for equation 
(85) are the pulsons. Pulsons are given by the sum over N velocity profiles determined by the 
Green's function g, as 

TV 

u{x,t) = J2Pa{t)g{x,qa{t)). (86) 

a=l 

Again for (85), the singular momentum ansatz (61) results in a finite-dimensional invariant mani- 
fold of solutions, whose dynamics is canonically Hamiltonian. The Hamiltonian for the canonical 
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dynamics of the 2A'" parameters Pa{t) and qa{t) in the "pulson" solutions (86) of equation (85) is 

1 ^ 

P<^Pb9{Qa,qb) ■ (87) 

a,b=l 

Again for the pulsons, the canonical equations for the invariant manifold of singular momentum 
solutions provide a phase-space description of geodesic motion, this time with respect to the co- 
metric given by the Green's function g. Mathematical analysis and numerical results for the 
dynamics of these pulson solutions are given in [26]. These results describe how the collisions of 
pulsons (86) depend upon their shape. 



Compactons in the l/a^ — )■ limit of CH. 

the CH equation (60) becomes the KdV equation, 
the Hunter-Zheng equation [36] 



As mentioned earlier, in the limit that — > 0, 
In the opposite hmit that l/a^ — > CH becomes 



{ut + uu^)^^ = ^{ul)x (Hunter-Zheng) 

This equation has "compacton" solutions, whose collision dynamics was studied numerically and 
put into the present context in [26]. The corresponding Green's function satisfies —dlg{x) — 25{x), 
so it has the triangTilar shape, g(x) = 1 — \x\ for \x\ < 1, and vanishes otherwise, for |a;| > 1. That 
is, the Green's function in this case has compact support; hence, the name "compactons" for these 
pulson solutions, which as a limit of the integrable CH equations are true solitons, solvable by 1ST. 



Pulson solutions of the dispersionless 6-equation. Holm and Staley [35] give the pulson 
solutions of the traveling wave problem and their elastic collision properties for the dispersionless 
6- equation, 

rrit + urrix + b mUx = , with u = g *m , (88) 

with any (symmetric) Green's function g and for any value of the parameter b. Numerically, 
pulsons and peakons are both found to be stable for 6 > 1, [35]. The reduction to noncanonical 
Hamiltonian dynamics for the invariant manifold of singular momentum solutions (61) of the other 
integrable case b — 3 with peakon Green's function g{x,y) = e~l^~^l/"^ is found in [19]. 



8 Euler-Poincare theory in higher dimensions 

Generalizing the peeikon solutions of the CH equation to higher dimensions. In [35], 
weakly nonlinear analysis and the assumption of columnar motion in the variational principle for 
Euler's equations were found to produce the two-dimensional generalization of the dispersionless 
CH equation (62). This generalization is the Euler-Poincare (EP) equation [34] for the Lagrangian 
consisting of the kinetic energy, 

dxdy , (89) 

in which the fluid velocity u is a two-dimensional vector. Evolution generated by kinetic energy 
in Hamilton's principle results in geodesic motion, with respect to the velocity norm ||u||, which 
is provided by the kinetic energy Lagrangian. For ideal incompressible fluids governed by Euler's 



]u]^ -I- Q;^(div u)' 
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equations, the importance of geodesic flow was recognized by Arnold [4] for the norm of the 
fluid velocity. The EP equation generated by any choice of kinetic energy norm without impos- 
ing incompressibility is called "EPDiff," for "Euler-Poincare equation for geodesic motion on the 
diffeomorphisms." EPDiff is given by [34] 



^— + u ■ V j m + Vu^ • m + m(div u) = , 



(90) 



with momentum density m — S£/5u , where i — |||u|p is given by the kinetic energy, which deflnes 
a norm in the fluid velocity ||u||, yet to be determined. By design, this equation has no contribution 
from either potential energy, or pressure. It conserves the velocity norm ||u|| given by the kinetic 
energy. Its evolution describes geodesic motion on the diffeomorphisms with respect to this norm 
[34]. An alternative way of writing the EPDiff equation (90) in either two, or three dimensions is. 



d_ 
Ft 



m — u X curl m + V(u • m) + m(div u) = . 



(91) 



This form of EPDiff involves all three differential operators, curl, gradient and divergence. For the 
kinetic energy Lagrangian i given in (89), which is a norm for irrotational flow (with curlu = 0), 
we have the EPDiff equation (90) with momentum m = 5£/5u = u — Q;^V(divu). 
EPDiff (90) may also be written intrinsically as 



dt 5u 



— — ad; 



(92) 



where ad* is the dual of the ad-operation (commutator) for vector fields. See [5, 43] for additional 
discussions of the beautiful geometry underlying this equation. 

Reduction to the dispersionless CH equation in ID. In one dimension, the EPDiff equation 
(90-92) with Lagrangian i given in (89) simplifies to the dispersionless CH equation (62). The 
dispersionless limit of the CH equation appears, because we have ignored potential energy and 
pressure. 



Strengthening the kinetic energy norm to allow for circulation. The kinetic energy 
Lagrangian (89) is a norm for irrotational fiow, with curlu = 0. However, inclusion of rotational 
fiow requires the kinetic energy norm to be strengthened to the norm of the velocity, defined 
as 



up-Fa^(divu)^-Fa^ (curl u) 

1 



dxdy 



lup -Fa^lVup 



dxdy = -||u||^i . 



(93) 



Here we assume boundary conditions that give no contributions upon integrating by parts. The 
corresponding EPDiff equation is (90) with m = S£/Su = u — Au . This expression involves 
inversion of the familiar Helmholtz operator in the (nonlocal) relation between fluid velocity and 
momentum density. The norm ||u||^i for the kinetic energy (93) also arises in three dimensions 
for turbulence modeling based on Lagrangian averaging and using Taylor's hypothesis that the 
turbulent fluctuations are "frozen" into the Lagrangian mean flow [22] . 
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Generalizing the CH peakon solutions to n dimensions. Building on the peakon solutions 
(59) for the CH equation and the pulsons (86) for its generalization to other traveling-wave shapes 
in [26], Holm and Staley [35] introduced the following measure- valued singular momentum solution 
ansatz for the ra— dimensional solutions of the EPDiff equation (90):^ 

^ r 

m(x,t) = ^ / P"(s,t)(5(x-Q"(s,t))rfs. (94) 

a=l ^ 

These singular momentum solutions, called "diffeons," are vector density functions supported in 
on a set of surfaces (or curves) of codimension {n — k) for s G M'^ with k < n. They may, 
for example, be supported on sets of points (vector peakons, k = 0), one-dimensional filaments 
(strings, k = 1), or two-dimensional surfaces (sheets, /c = 2) in three dimensions. 

Figure 3 shows the results for the EPDiff equation when a straight peakon segment of finite 
length is created initially moving rightward (East). Because of propagation along the segment in 
adjusting to the condition of zero speed at its ends and finite speed in its interior, the initially 
straight segment expands outward as it propagates and curves into a peakon "bubble." 




Figure 3: A peakon segment of finite length is initially moving rightward (East). Because its speed 
vanishes at its ends and it has fully two-dimensional spatial dependence, it expands into a peakon 
"bubble" as it propagates. (The colors indicate speed: red is highest, yellow is less, blue low, grey 
zero.) 

Figure 4 shows an initially straight segment whose velocity distribution is exponential in the 
transverse direction, but is wider than a for the peakon solution. This initial velocity distribution 
evolves under EPDiff to separate into a train of curved peakon "bubbles," each of width a. This 
example illustrates the emergent property of the peakon solutions in two dimensions. This phe- 
nomenon is observed in nature, for example, as trains of internal wave fronts in the south China 
Sea [42]. 

^These solutions represent smooth embeddings Emb(M'^,]R") with k < n. In contrast, the similar expression 
(46) for the soliton solutions represent smooth functions M -> M. 
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Figure 4: An initially straight segment of velocity distribution whose exponential profile is wider 
than the width a for the peakon solution will break up into a train of curved peakon "bubbles," 
each of width a. This example illustrates the emergent property of the peakon solutions in two 
dimensions. 



Substitution of the singular momentum solution ansatz (94) into the EPDiff equation (90) im- 
plies the following integro-partial-differential equations (IPDEs) for the evolution of the parameters 
{P} and {Q}, 



d_ 
di 

d_ 
di 



Q"(^'^) = /PV,i)G(Q"(s,t)-Q^(s',t))cis' 
b=i 

N 

L 1 ^ 



(95) 



d 



G{Q''{s,t)-Cl''{s',t))ds' . 



Importantly for the interpretation of these solutions, the coordinates s G M'^ turn out to be 
Lagrangian coordinates. The velocity field corresponding to the momentum solution ansatz (94) 
is given by 



u(x, t) = G * m 



N 

J2 / P\s',t)G{^-Q\s',t))ds'. 
b=i 



(96) 
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Figure 5: A single collision is shown involving reconnection as the faster peakon segment initially 
moving Southeast along the diagonal expands, curves and obliquely overtakes the slower peakon 
segment initially moving rightward (East). This reconnection illustrates one of the collision rules 
for the strongly two-dimensional EPDiff flow. 



for u G M". When evaluated along the curve x = Q"(s, t), this velocity satisfies, 

N 

b=i 

dQ%s,t) 
dt ■ 



(97) 



Consequently, the lower-dimensional support sets defined on x = Q"'{s,t) and parameterized by 
coordinates s G move with the fluid velocity. This means the s G M'^ are Lagrangian coordi- 
nates. Moreover, equations (95) for the evolution of these support sets are canonical Hamiltonian 
equations, 

^Q°(^.*) = ^. |P°(^.*) = -S^- (98) 



dt^ ' ' SP" dt 
The corresponding Hamiltonian function Hj\f : (M" x M")^ 



K IS, 



H 



N 



N 



•^•^ a,b=l 

G{Ct{s,t),Q\s',t))ds ds' . 



(99) 



This is the Hamiltonian for geodesic motion on the cotangent bundle of a set of curves Q"(s,t) 
with respect to the metric given by G. This dynamics was investigated numerically in [35] to 
which we refer for more details of the solution properties. One important result found numerically 
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in [35] is that only codimension-one singular momentum solutions appear to be stable under the 
evolution of the EPDiff equation. Thus, we have 

Stability for codimension-one: the singular momentum solutions of EPDiff' are stable, as points 
on the line (peakons), as curves in the plane (filaments, or wave fronts), or as surfaces in space 
(sheets). 

Proving this stability result analytically remains an outstanding problem. The stability of 
peakons on the real line is proven in [17]. 

Reconnections in oblique overtaking collisions of peakon wave fronts. Figures 5 shows 
results of oblique wave front collisions producing reconnections for the EPDiff equation in two 
dimensions. Figure 5 shows a single oblique overtaking collision, as a faster expanding peakon 
wave front overtakes a slower one and reconnects with it at the collision point. 

The pecikon solution ansatz is a momentum map. As shown in [33] , the solution expressions 
(61) in one dimension and (94) in higher dimensions may be interpreted as equivariant momentum 
maps, from the cotangent bundle of the smooth embeddings of lower dimensional sets M'* C M" , to 
the dual of the Lie algebra of vector fields defined on these sets. (Momentum maps for Hamiltonian 
dynamics are reviewed in [43], for example.) The result that the singular solution ansatz (94) is 
a momentum map helps to organize the theory, to explain previous results and to suggest new 
avenues of exploration. This geometric feature underlies the remarkable reduction properties of 
the CH and EPDiff equations, and explains why they must be Lie-Poisson Hamiltonian equations. 
This is because of the general fact that equivariant momentum maps are Poisson maps. This 
geometric feature also underlies the singular momentum solution (94) and its associated velocity 
(96) which generalize the peakon solutions, both to higher dimensions and to arbitrary kinetic 
energy metrics. As we saw in Section 3, the soliton solution (53) is also a momentum map. This 
soliton momentum map may be expected to apply in the action-angle representation of the solution 
of any integrable Hamiltonian PDE. Its further properties will be studied in detail elsewhere. 

9 Three open problems 

(1.) Throughout this discussion the solutions u{x, t) were confined to be functions in the Schwartz 
class, a; > 0. The situation when the condition m{x, 0) + w > on the initial data does not hold is 
more complicated and requires separate analysis [39, 10, 12]. In general, it leads to wave-breaking 
[12]. An attempt at developing the inverse scattering theory for this case has been made by Kaup 
[39], who suggested applying the inverse scattering approach separately in each interval where 
m{x, t) + uj is of the same sign. The problem however is how to join solutions that are valid in 
different intervals. 

(2.) The peakon solution (68) was obtained from the soliton solution under the assumption 
m{x, 0) + a; > 0. Thus, all pk are of the same sign, since all the eigenvalues A„ in this case are 
positive. However, one can formally use the same solution with eigenvalues of various signs to 
model Pk of various signs (mixture of peakons and 'anti-peakons') and thus to study peakon - 
anti-peakon interactions, see e.g. [51]. The result is that the multi-peakon interaction (including 
anti-peakons) in general decomposes into a sequence of pairwise collisions [26]. The collision of a 
single peakon-anti-peakon pair was studied already in [9] . When the eigenvalues are of mixed signs, 
the Hankel determinants in the denominator of (68) may develop singularities for finite values of 
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t. This 'peakon-breaking' phenomenon is apparently the analog of the wave-breaking mentioned 
earlier when — ^ 0. This needs to be investigated further. 

(3.) Stability for EPDifT singular momentum solutions, that is, proving their stability analyti- 
cally, remains an outstanding problem. 
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